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We investigate the temperature dependence of the coefficients in the Landau-Lifshitz-Gilbert 
equation in ferromagnetic GaMnAs by employing the Zener model. We first calculate the hole spin 
relaxation time based on the microscopic kinetic equation. We find that the hole spin relaxation 
time is typically several tens femtoseconds and can present a nonmonotonic temperature dependence 
due to the variation of the interband spin mixing, influenced by the temperature related Zeeman 
splitting. With the hole spin relaxation time, we are able to calculate the coefficients in the Landau- 
Lifshitz-Gilbert equation, such as the Gilbert damping, nonadiabatic spin torque, spin stiffness and 
vertical spin stiffness coefficients. We find that the nonadiabatic spin torque coefficient j3 is around 
0.1 ~ 0.3 at low temperature, which is consistent with the experiment [Adam et al, Phys. Rev. B 
80, 193204 (2009)]. As the temperature increases, /3 monotonically increases and can exceed one 
in the vicinity of the Curie temperature. In the low temperature regime with /3 < 1, the Gilbert 
damping coefficient a increases with temperature, showing good agreement with the experiments 
[Sinova et al, Phys. Rev. B 69, 085209 (2004); Khazen et al, ibid. 78, 195210 (2008)]. Furthermore, 
we predict that a decreases with increasing temperature once ft > 1 near the Curie temperature. 
We also find that the spin stiffness decreases with increasing temperature, especially near the Curie 
temperature due to the modification of the finite ft. Similar to the Gilbert damping, the vertical 
spin stiffness coefficient is also found to be nonmonotonically dependent on the temperature. 

PACS numbers: 72.25.Rb, 75.50.Pp, 72.25.Dc, 75.30.Gw 



I. INTRODUCTION 

The ferromagnetic semiconductor, GaMnAs, has been 
proposed to be a promising candidate to realize all- 
semiconductor spintronic devices r* 2 - where the existence 
of the ferromagnetic phase in the heavily doped sample 
sustains seamless spin injection and detection in normal 
non-magnetic semiconductors. 3,4 One important issue for 
such applications lies in the efficiency of the manipula- 
tion of the macroscopic magnetization, which relies on 
properties of the magnetization dynamics. Theoretically, 
the magnetization dynamics can be described by the ex- 
tended Landau-Lifshitz-Gilbert (LLG) equation^ - — 

ri = —711 x H c ff + an x h — (1 — /3nx)(v s • V)n 

-inx(4-i>x)V 2 n, (1) 

with n and Md standing for the direction and magni- 
tude of the magnetization, respectively. H e g is the ef- 
fective magnetic field and/or the external field. The sec- 
ond term on the right hand side of the equation is the 
Gilbert damping torque with a denoting the damping 
coefficient^'^ The third one describes the spin-transfer 
torque induced by the spin current v s 2& As reported, 
the out-of-plane contribution of the spin-transfer torque, 
measured by the nonadiabatic torque coefficient /3, can 
significantly ease the domain wall motion.— £ In Eq. (fT|), 
the spin stiffness and vertical spin stiffness coefficients are 
evaluated by A ss and A^ s respectively, which are essen- 
tially important for the static structure of the magnetic 
domain walli 1 ^ Therefore, for a thorough understanding 



of properties of the magnetization dynamics, the exact 
values of the above coefficients are required. 

In the past decade, the Gilbert damping and nonadia- 
batic torque coefficients have been derived via many mi- 
croscopic approaches, such as the Blotzmann equation, 11 
diagrammatic calculation ) 12 ' 13 Fermi-surface breathing 
modeli 4 ^— and kinetic spin Bloch equations! 10 ' 17 Accord- 
ing to these works, the spin lifetime of the carriers was 
found to be critical to both a and /3. However, to the 
best of our knowledge, the microscopic calculation of the 
hole spin lifetime in ferromagnetic GaMnAs is still ab- 
sent in the literature, which prevents the determination 
of the values of a and j3 from the analytical formulas. Al- 
ternatively, Sinova et al^& identified the Gilbert damp- 
ing from the susceptibility diagram of the linear-response 
theory and calculated a as function of the quasiparticle 
lifetime and the hole density. Similar microscopic calcu- 
lation on j3 was later given by Garate et al.J£ In those 
works, the quasiparticle lifetime was also treated as a 
parameter instead of explicit calculation. Actually, the 
accurate calculation of the hole spin and/or quasiparti- 
cle lifetime in ferromagnetic GaMnAs is difficult due to 
the complex band structure of the valence bands. In the 
present work, we employ the microscopic kinetic equation 
to calculate the spin lifetime of the hole gas and then eval- 
uate a and f3 in ferromagnetic GaMnAs. For the velocity 
of the domain-wall motion due to the spin current, the 
ratio /3/a is an important parameter, which has attracted 
much attention . 12 ' 19 ' 20 Recently, a huge ratio (~ 100) in 
nanowire was predicted from the calculation of the scat- 
tering matrix by Hals et al.£Q By calculating a and /3, we 
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are able to supply detailed information of this interesting 
ratio in bulk material. Moreover, the peak-to-peak fer- 
romagnetic resonance measurement revealed pronounced 
temperature and sample preparation dependences of the 
Gilbert damping coefficient] 18 ' 21 ' 22 For example, in an- 
nealed samples, a can present an increase in the vicinity 
of the Curie temperature ; 18 ! 21 which has not been stud- 
ied theoretically in the literature. Here, we expect to 
uncover the underlying physics of these features. In ad- 
dition, the nonadiabatic torque coefficient j3 in GaMnAs 
has been experimentally determined from the domain- 
wall motion and quite different values were reported by 
different groups, from 0.01 to 0.36 ) 23 i 24 which need to be 
verified by the microscopic calculation also. Moreover, 
to the best of our knowledge, the temperature depen- 
dence of j3 has not been studied theoretically. We will 
also address this issue in the present work. 

In the literature, the spin stiffness in GaMnAs was 
studied by Konig et aL ) 25 ' 26 who found that A ss increases 
with hole density due to the stronger carrier-mediated 
interaction between magnetic ions, i.e., A ss = A f ^/(4m*) 
with Nh and m* being the density and effective mass 
of hole gas, separately. However, as shown in our pre- 
vious work, the stiffness should be modified as A ss ~ 
A r / 1 /[4m*(l + /3 2 )] in ferromagnetic GaMnAs with a finite 
/Jiifi As a result, A ss as well as the vertical spin stiffness 
A^ s = /3A SS may show a temperature dependence intro- 
duced by j3. This is also a goal of the present work. 

For a microscopic investigation of the hole dynamics, 
the valence band structure is required for the descrip- 
tion of the occupied carrier states. In the literature, 
the Zener model^l based on the mean-field theory has 
been widely used for itinerant holes in GaMnAs^Sr— 
where the valence bands split due to the mean-field p- 
d exchange interaction. In the present work, we utilize 
this model to calculate the band structure with the ef- 
fective Mn concentration from the experimental value of 
the low-temperature saturate magnetization in GaMnAs. 
The thermal effect on the band structure is introduced 
via the temperature dependence of the magnetization fol- 
lowing the Brillouin function. Then we obtain the hole 
spin relaxation time by numerically solving the micro- 
scopic kinetic equations with the relevant hole-impurity 
and holc-phonon scatterings. The carrier-carrier scatter- 
ing is neglected here by considering the strongly degen- 
erate distribution of the hole gas below the Curie tem- 
perature. We find that the hole spin relaxation time 
decreases/increases with increasing temperature in the 
small/large Zeeman splitting regime, which mainly re- 
sults from the variation of the interband spin mixing. 
Then we study the temperature dependence of the co- 
efficients in the LLG equation, i.e., a, /3, A ss and A^ s , 
by using the analytical formulas derived in our previous 
works. 10 ' 17 Specifically, we find that /3 increases with in- 
creasing temperature and can exceed one in the vicinity 
of the critical point, resulting in very interesting behav- 
iors of other coefficients. For example, a can present an 
interesting nonmonotonic temperature dependence with 



the crossover occurring at (3 ~ 1. Specifically, a increases 
with temperature in the low temperature regime, which is 
consistent with the experiments. Near the Curie temper- 
ature, an opposite temperature dependence of a is pre- 
dicted. Similar nonmonotonic behavior is also predicted 
in the temperature dependence of A^ s . Our results of /3 
and ^4 SS also show good agreement with the experiments. 

This work is organized as follows. In Sec. II, we setup 
our model and lay out the formulism. Then we show 
the band structure from the Zener model and the hole 
spin relaxation time from microscopic kinetic equations 
in Sec. III. The temperature dependence of the Gilbert 
damping, nonadiabatic spin torque, spin stiffness and 
vertical spin stiffness coefficients are also shown in this 
section. Finally, we summarize in Sec. IV. 

II. MODEL AND FORMULISM 

In the sp-d model, the Hamiltonian of hole gas in GaM- 
nAs is given by 3 -^ 

H = H p + H p d, (2) 

with Hp describing the itinerant holes. H p d is the sp-d 
exchange coupling. By assuming that the momentum k 
is still a good quantum number for itinerant hole states, 
one employs the Zener model and utilizes the k • p pertur- 
bation Hamiltonian to describe the valence band states. 
Specifically, we take the eight-band Kane Hamiltonian 
i/jf(k) fRef.l32T) in the present work. The sp-d exchange 
interaction reads 

Hp d = -j^yJ2 E J ex m 'S* • <mk|i|m'k)4 kCWk , 

I mm'k 

(3) 

with No and V standing for the density of cation sites 
and the volume, respectively. The cation density Nq = 
2.22 x 10 22 cm -3 . The eight-band spin operator can be 
written as J = (|<r) © J 3 / 2 © J1/2, where |er, J 3 / 2 and 
J 1/2 represent the total angular momentum operators of 
the conduction band, Tg valence band and Tf valence 
band, respectively. J^™ 1 ' stands for the matrix element 
of the exchange coupling, with {m} and {m'} being the 
basis defined as the eigenstates of the angular momen- 
tum operators J . The summation of is through all 
localized Mn spins S; (at r;). 

Then we treat the localized Mn spin in a mean-field 
approximation and obtain 

H pd = -x eS {S) • I J2 ^ ro Vk|i|m'k)4 kCm , k J , 

(4) 

where (S) represents the average spin polarization of 
Mn atoms with uncompensated doping density Nmh = 
XeflNo. Obviously, H p d can be reduced into three blocks 
as J, i.e., fl"^ m '(k) = A mm n- (mk\J\m'k) with the Zee- 
man splitting of the m-band A mm = -x eS S d J^ m ^jm- 
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Here, n is the direction of (S). For a manganese ion, the 
total spin Sd = 5/2. The temperature-dependent spon- 
taneous magnetization M(T) can be obtained from the 
following equation of the Brillouin function^ 



BaM = 



S d + 1 T 
35 d T L 



(5) 



MviTmT' with Tc being the Curie 



2S £j ±l coth( 2S^l y) 



2S d 



where y 

temperature. Here, Bs d {y) 
i, coth(^j/). 

The Schrddinger equation of the single particle Hamil- 
tonian is then written as 



[H K (k) + H pd (k)]\f, 7 k) = E^k). 



(6) 



One obtains the band structure and wave functions from 
the diagonalization scheme. In the presence of a finite 
Zeeman splitting, the structure of the valence bands devi- 
ates from the parabolic dispersion and becomes strongly 
anisotropic as we will show in the next section. Moreover, 
the valence bands at Fermi surface are well separated in 
ferromagnetic GaMnAs because of the high hole density 
(> 10 20 cm~ 3 ) and Zeeman splitting, suggesting that the 
Fermi golden rule can be used to calculate the lifetime of 
the quasiparticle states. For example, the contribution of 
the hole-impurity scattering on the /^th-band state with 
energy e can be expressed by 



x <5(6 Mk - 6, q )(7 2 _ q |( M k|^q)| 2 /(^k)[l - /(e*q)], (7) 

where -D M (e) stands for the density of states of the /xth 
band, /(e^k) satisfies the Fermi distribution in the equi- 
librium state. The hole-impurity scattering matrix ele- 
ment J7 2 = Z 2 e 4 /[K (q 2 + k 2 )] 2 with Z = 1. k and 
k denote the static dielectric constant and the screening 
constant under the random-phase approximation^ re- 
spectively. Similar expression can also be obtained for 
the hole-phonon scattering. 

However, it is very complicated to carry out the multi- 
fold integrals in Eq. ([7]) numerically for an anisotropic 
dispersion. Also the lifetime of the quasiparticle is not 
equivalent to the spin lifetime of the whole system, which 
is required to calculate the LLG coefficients according 
to our previous work . 10 ' 17 Therefore, we extend our ki- 
netic spin Bloch equation approach 3 - to the current sys- 
tem to study the relaxation of the total spin polarization 
as follows. By taking into account the finite separation 
between different bands, one neglects the interband co- 
herence and focuses on the carrier dynamics of the non- 
equilibrium population. The microscopic kinetic equa- 
tion is then given by 



9*n M ,k = <9*?V,k| 



d t n 



/i,k | 



I hp 



(8) 



with n Ml k being the carrier occupation factor at the /ith 
band with momentum k. The first and second terms on 



the right hand side stand for the hole-impurity and hole- 
phonon scatterings, respectively. Their expressions can 
be written as 



d t n 



-27m, E/ k ! -k'('Vk - n vk ,)\(fik\vk')\ 2 



- k' 



^{E^ — Ey k >), 



and 



d t n 



I hp 



-2tt 



E 

\,±,v,k> 



(9) 



'8{E uk l — E^k ± UJ\,q) 



i±i. The details of 



with N x,i = [exp(w A , q /fcsT)-l] 
the hole-phonon scattering elements |Mq| 2 can be found 



in Refs. l36l - l38l From an initial condition with a small 
non-equilibrium distribution, the temporal evolution of 
the hole spin polarization is carried out by 



J(t) 



-E 

/j.,k 



(fik\J\iik)n^(t), 



(11) 



from the numerical solution of Eq. ((HJ. The hole spin 
relaxation time can be extracted from the exponential 
fitting of J with respect to time. One further calculates 
the concerned coefficients such as a, /3, A ss and Al s . 



III. NUMERICAL RESULTS 

In the Zener model, the sp-d exchange interaction con- 
stants J^ m are important parameters for the band struc- 
ture. In the experimental works, the p-d exchange cou- 
pling constant Jpp was reported to vary from —1 eV 
to 2.5 eV, depending on the doping density^2r— In 
ferromagnetic samples, JJ^ is believed to be negative, 
which was demonstrated by theoretical estimation Jpp 
—0.3 eV (Ref.H2). In our calculation, the antiferromag- 
netic p-d interaction JPg is chosen to be —0.5 eV or 
— 1.0 eV. The ferromagnetic s-d exchange coupling con- 
stant is taken to be = 0.2 eV (Ref.l3lT). 

Another important quantity for determining the Zee- 
man splitting is the macroscopic magnetization or the 
effective concentration of the Mn atoms. As deduced 
from the low-temperature saturate magnetization, only 
around 50% Mn atoms can contribute to the ferromag- 
netic magnetization, which has been recognized as the in- 
fluence of the compensation effect due to the deep donors 
(e.g., As antisites) or the formation of sixfold-coordinated 
centers defect Mn 6As (RefJUj). As only the uncompen- 
sated Mn atoms can supply holes and contribute to the 
ferromagnetic magnetization^ one can also estimate the 
total hole density from the saturate magnetization. 45 
However, the density of the itinerant hole can be smaller 
than the effective Mn concentration because of the local- 
ized effect in such disordered material. It was reported 



4 





T c 
(K) 


M, 
(emu-cm -3 ) 


NMn 

(10 20 cm" 3 ) 




130 


38 


8 


B a 


157 


47 


10 


C b 


114 


33 


6.9 


D c 


110 






E d 


139 


53.5 


11.5 



momentum space. 



Rcf. 21 



Rei.m 
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TABLE I: The parameters obtained from the experiments 
for different samples: A: Gao.93Mno.07As/Gao.902Ino.09sAs; B: 
Gao.93Mno.07As/GaAs; C: Gao.93Mno.o7As/Gai_j,In y As; D: 
Gao.92Mno.08As; E: Gao.896Mno.104Aso.93Po. 07- M s stands for 
the saturate magnetization at zero temperature M(0). 



that the hole density is only 15-30% of the total concen- 
tration of the Mn atoms.— 

In our calculation, the magnetization lies along the 
principle axis chosen as [001]-direction. 3 A The conven- 
tional parameters are mainly taken from those of GaAs 
in Refs.H^ and H3- Other sample-dependent parame- 
ters such as the Curie temperature and effective Mn 
concentration are picked up from the experimental 
works For sample A, B and E (C), only the 
saturate magnetization at 4 (104) K was given in the 
references. Nevertheless, one can extrapolate the zero 
temperature magnetization M s from Eq. ([5]). The effec- 
tive Mn concentrations listed in Table|T] are derived from 
A^Mn = M s / (gfisSd)- It is clear to see that all of these ef- 
fective Mn concentrations are much smaller than the dop- 
ing density (> 1.5 x 10 21 cm -3 ) due to the compensation 
effect as discussed above. Since the saturate magneti- 
zation of sample D is unavailable, we treat the effective 
Mn concentration as a parameter in this case. More- 
over, since the exact values of the itinerant hole densities 
are unclear in such strongly disordered samples, we treat 
them as parameters. Two typical values are chosen in 
our numerical calculation, i.e., Nh — 3 x 10 20 cm -3 and 
5 x 10 20 cm -3 . The effective impurity density is taken to 
be equal to the itinerant hole density. 

For numerical calculation of the hole spin dynamics, 
the momentum space is partitioned into blocks. Com- 
pared to the isotropic parabolic dispersion, the band 
structure in ferromagnetic GaMnAs is much more com- 
plex as we mentioned above [referred to Figs.QJb) and 
5]. Therefore, we need to extend the partition scheme 
used in isotropic parabolic dispersion^ into anisotropic 
case. In our scheme, the radial partition is still carried 
out with respect to the equal-energy shells, while the an- 
gular partition is done by following Ref.|48|. In contrast 
to the isotropic case, the number of states in one block is 
generally different from that in another block even both 
of them are on the same equal-energy shell. We calculate 
the number of states of each block from its volume in 
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FIG. 1: (Color online) (a) Zeeman energy as function of tem- 
perature, (b) The valence band structure with A pp = 45 meV. 
The blue dashed curve illustrates the Fermi level for the hole 
density Nh = 3 x 10 20 cm -3 , while the green dotted one gives 
Nh = 5 x 10 20 cm -3 . The density of states as function of 
energy at (c) T/T c = 0.1 and (d) T/T c = 0.99 for the uncom- 
pensated Mn density NtAn = 8 x 10 21 cm~ 3 . In (d), the blue 
dashed curve stands for the upper heavy hole band from the 
spherical approximation and the corresponding DOS from the 
analytical formula {-\/2E[\/m* / '(27rft)] 3 ) is given as the green 
dotted curve. Here, JS = —0.5 eV. 



A. Density of states 

By solving Eq. ([5]) , one obtains the magnetization at 
finite temperature M(T) and the corresponding Zeeman 
energy A pp . In Fig.QJa), the Zeeman splitting from 
-^ex = — 0.5 eV is plotted as function of the temperature. 
It is seen that the Zeeman energy is tens of milli-electron 
volts at low temperature and decreases sharply near the 
Curie temperature due to the decrease of the magnetiza- 
tion. To show the anisotropic nonparabolic feature of the 
band structure in the presence of the Zeeman splitting, 
we illustrate the valence bands along [001]- and [In- 
directions in Fig.QJb), which are obtained from Eq. (|6]) 
at T/T c = 0.1 for A^ Mn = 8 x 10 20 cm~ 3 . In this case, the 
Zeeman splitting A pp = 45 meV. The Fermi levels for the 
hole densities Nh — 3 x 10 20 cm -3 and 5 x 10 20 cm -3 are 
shown as blue dashed and green dotted curves, respec- 
tively. As one can see that all of the four upper bands 
can be occupied and the effective mass approximation 
obviously breaks down. 

By integrating over the volume of each equal-energy 
shell, one obtains the density of states (DOS) of each 
band as function of energy in Fig.QJc) and (d). Here the 
energy is defined in the hole picture so that the sign of 
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the energy is opposite to that in Fig.[ljb). It is seen that 
the DOS of the upper two bands are much larger than 
those of the other bands, regardless of the magnitude of 
the Zeeman splitting. For T/T c — 0.99, the systems ap- 
proaches the paramagnetic phase and the nonparabolic 
effect is still clearly seen from the DOS in Fig.QJd), es- 
pecially in the high energy regime. Moreover, the pro- 
nounced discrepancy of the DOS for the two heavy hole 
bands suggests the finite splitting between them. We 
find that these features are closely connected with the 
anisotropy of the valence bands, corresponding to the 
Luttinger parameters 72 ^ 73 in GaAs4^ In our calcu- 
lation, we take 71 = 6.85, 72 = 2.1 and 73 = 2.9 from 
Ref.l47l. As a comparison, we apply a spherical approx- 
imation (71 = 6.85 and 72 = 73 = 7 = 2.5) and find 
that the two heavy hole bands become approximately 
degenerate.— The DOS of the upper heavy hole band 
is shown as the blue dashed curve in Fig.QJd), where 
we also plot the corresponding DOS from the analyti- 
cal expression, i.e., V2E[\^m* / (2nh)} 3 , as the green dot- 
ted curve. Here, we use the heavy-hole effective mass 
m* = m 1 (7! — 27) with mo denoting the free electron 
mass. The perfect agreement between the analytical and 
our numerical results under the spherical approximation 
suggests the good precision of our numerical scheme. 
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FIG. 2: (Color online) The equilibrium hole spin polarization 
as function of Zeeman splitting for sample A and B. Here, 
JEi = -0.5 eV. 



B. Hole spin relaxation 

In this part, we investigate the hole spin dynamics by 
numerically solving the microscopic kinetic equation, i.e., 
Eq. ([5]) . By taking into account the equilibrium hole spin 
polarization, we fit the temporal evolution of the total 
spin polarization along [001]-direction by 

JS) = J% + J' z e- t/T % (12) 

where J® and J' z correspond to the equilibrium and non- 
equilibrium spin polarizations, respectively. r s is the hole 
spin relaxation time. 



In all the cases of the present work, the equilibrium 
hole spin polarization for a fixed hole density is found 
to be approximately linearly dependent on the Zeeman 
splitting. In Fig. [3J 3® in samples A and B (similar be- 
havior for others) are plotted as function of Zeeman split- 
ting, where the exchange coupling constant is taken 
to be —0.5 eV. One notices that the average spin polar- 
ization becomes smaller with the increase of the hole den- 
sity, reflecting the large interband mixing for the states 
in the high energy regime. 
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FIG. 3: (Color online) (a) Spin relaxation time as function of 
temperature with = —0.5 eV for sample A and B. (b) Spin 
relaxation time as function of temperature and Zeeman split- 
ting obtained from the calculation with J£f = — 1 eV for sam- 
ple B. The inset at the left (right) upper corner illustrates the 
band structure from [001]-direction to [11 Indirection [refer to 
Fig-DIb)] for T/T c = 0.4 (0.99) and A pp = 105 (16.7) meV. 
The Fermi levels of N h = 3 x 10 20 cm" 3 and 5 x 10 20 cm" 3 
are shown as the blue dashed and green dotted curves in the 
insets, separately. 



The temperature dependence of the hole spin relax- 
ation time in samples A and B with JJ^ = —0.5 eV is 
shown in Fig.^a), where the spin relaxation time mono- 
tonically decreases with increasing temperature. This 
feature can be understood from the enhancement of 
the interband mixing as the Zeeman splitting decreases 
(shown below) £H To gain a complete picture of the role of 
the Zeeman splitting on the hole spin relaxation in fer- 
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romagnetic GaMnAs, we also carry out the calculation 
with the exchange constant J pp — — 1 eVi 31 i 39 Very in- 
terestingly, one finds that the hole spin relaxation time at 
low temperature increases with increasing temperature, 
resulting in a nonmonotonic temperature dependence of 
the hole spin relaxation time in sample B. The results 
in this case are shown as solid curves in Fig.[3Jb), where 
we also plot the Zeeman splitting dependence of the hole 
spin relaxation time as dashed curves. It is seen that 
the hole spin relaxation time for the hole density Nh — 
3 x 10 20 cm -3 first increases with increasing temperature 
(alternatively speaking, decreasing Zeeman splitting) and 
starts to decrease at around 0.8 T c where the Zeeman 
splitting A pp — 70 meV. To understand this feature, we 
show the typical band structure in the increase (decrease) 
regime of the hole relaxation time at T/T c — 0.4 (0.99), 
corresponding to A pp = 105 (16.7) meV, in the inset at 
the left (right) upper corner. The Fermi levels of the 
hole density 3 x 10 20 cm~ 3 are labeled by blue dashed 
curves. One finds that the carrier occupations in the 
increase and decrease regimes are quite different. Specif- 
ically, all of the four upper bands are occupied in the 
decrease regime while only three valence bands are rele- 
vant in the increase regime. 

One may naturally expect that the increase regime 
originates from the contribution of the fourth band via 
the inclusion of the additional scattering channels or the 
modification of the screening. However, we rule out this 
possibility through the computation with the fourth band 
artificially excluded, where the results are qualitatively 
the same as those in Fig.(3{b). Moreover, the variations 
of the screening and the equilibrium distribution at fi- 
nite temperature are also demonstrated to be irrelevant 
to the present nonmonotonic dependence by our calcula- 
tion (not shown here). Therefore, the interesting feature 
has to be attributed to the variations of the band dis- 
tortion and spin mixing due to the exchange interaction. 
This is supported by our numerical calculation, where 
the nonmonotonic behavior disappears once the effect of 
the interband mixing is excluded by removing the wave- 
function overlaps |(^k|^k')| 2 in Eqs. (HJ) and (fTO]) (not 
shown here). 

For a qualitative understanding of the nonmonotonic 
temperature dependence of the hole spin relaxation time, 
we plot the Fermi surface in the k x -k z (k y = 0) and k x - 
k y (k z = 0) planes at N h = 3 x 10 20 cm" 3 in FigJU We 
choose typical Zeeman splittings in the increase regime 
(A pp = 105 meV), the decrease regime {A pp = 16.7 meV) 
and also the crossover regime (A pp = 70 meV). One no- 
tices that the Fermi surfaces in Fig. Ufa) and (d) are com- 
posed of three closed curves, meaning that only three 
bands are occupied for A pp = 105 meV [also see the in- 
set of Fig.|3{b)]. For the others with smaller Zeeman 
splittings, all of the four upper bands are occupied. The 
spin expectation of each state at Fermi surface is repre- 
sented by the color coding. Note that the spin expecta- 
tion of the innermost band for A pp — 70 meV is close to 
— 1.5 [see Fig.QJb) and (e)], suggesting that this band is 
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FIG. 4: (Color online) The Fermi surface in the k x -k z (k y — 0) 
and k x -k y (k z = 0) planes with A pp =105 meV (a,d), 70 meV 
(b,e) and 16.7 meV (c,f). The color coding represents the 
spin expectation of each state, £ = {^{Jzlfi) . Here, Nh = 
3 x 10 2 ° cm -3 . 



the spin-down heavy hole band and the mixing of other 
spin components in this band is marginal. Therefore, 
the spin-flip scattering related to this band is weak and 
can not result in the increase of the hole spin relaxation 
time mentioned above. By comparing the results with 
A pp — 105 meV and 70 meV, one notices that the spin 
expectation of the Fermi surface of the outermost band 
is insensitive to the Zeeman splitting. Therefore, this 
band can not be the reason of the increase regime also. 
Moreover, for the second and third bands in Fig.^Ja) 
and (b), from the comparable color coding between the 
two figures in this regime [also see Fig.HJd) and (e) with 
k z = 0], one finds that the spin expectation for the 
states with small k z is also insensitive to the Zeeman 
splitting. However, for the states with large k z , the spin 
expectation of the spin-down states (£ < 0) approaches 
a large magnitude (—1.5) with decreasing Zeeman split- 
ting, suggesting the decrease of the mixing from the spin- 
up states. As a result, the interband spin-flip scattering 
from/to these states becomes weak and the hole spin re- 
laxation time increases. In the decrease regime of the 
hole spin relaxation time, Fig.Htc) and (f) show that the 
two outer/inner bands approach each other, leading to a 
strong and anisotropic spin mixing. Therefore, the spin- 
flip scattering becomes more efficient in this regime and 
the spin relaxation time decreases. One may suppose 
that the nonmonotonic temperature dependence of the 
hole spin relaxation time can also arise from the varia- 
tion of the shape of the Fermi surface, according to Fig. 2) 
However, this variation itself is not the key of the non- 
monotonic behavior, because the calculation with this 
effect but without band mixing can not recover the non- 
monotonic feature as mentioned in the previous para- 
graph. For the hole density Nh = 5 x 10 20 cm~ 3 , the 
structures of the Fermi surface at A pp = 105 meV are 
similar to those in Fig.[4|b) and (e) . This explains the ab- 
sence of the increase regime for this density in Fig.[3]Jb). 
Moreover, we should point out that the increase regime 
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of the hole spin relaxation time in sample A for Jpp = 
— 1 eV is much narrower than that in sample B. The 
reason lies in the fact of lower effective Mn density in 
sample A, leading to the smaller maximal Zeeman split- 
ting ~ 90 meV, only slightly larger than the crossover 
value 70 meV at Nh = 3 X 10 20 cm" 3 . 

As a summary of this part, we find different temper- 
ature dependences of the hole spin relaxation time due 
to the different values of effective Mn concentration, hole 
density and exchange coupling constant J£g. In the case 
with large coupling constant and high effective Mn con- 
centration, the interband spin mixing can result in a non- 
monotonic temperature dependence of the hole spin re- 
laxation time. Our results suggest a possible way to esti- 
mate the exchange coupling constant with the knowledge 
of itinerant hole density, i.e., by measuring the temper- 
ature dependence of the hole spin relaxation time. Al- 
ternatively, the discrepancy between the hole relaxation 
time from different hole densities in Fig. [3][b) suggests 
that one can also estimate the itinerant hole density if 
the exchange coupling constant has been measured from 
other methods. 



Gilbert damping and non-adiabatic torque 
coefficients 



Facilitated with the knowledge of the hole spin re- 
laxation time, we can calculate the coefficients in the 
LLG equation. According to our previous works ; 10 ! 17 
the Gilbert damping and nonadiabatic spin torque co- 
efficients can be expressed as 



and 



a = J h /[N Mn \(S)\(/3 + l/l3)} 



(13) 



(14) 



respectively. In Eq. (TTBI . Jh represents the total equi- 
librium spin polarization of the itinerant hole gas, i.e., 
Jh = N h J° with J° being the one defined in Eq. (JT2J in 
our study. The average spin polarization of a single Mn 
ion is given by |(S)| = S d M(T)/M(0). 

In Fig.EJa), (c) and (e), the nonadiabatic spin torque 
coefficients (3 in sample A-C are plotted as function of 
temperature with Jpp = —0.5 eV and —1.0 eV. Our re- 
sults in sample C show good agreement with the experi- 
mental data (plotted as the brown square) in Fig.[5je). 23 
At low temperature, the value of j3 is around 0.1~0.3, 
which is also comparable with the previous theoretical 
calculation^ Very interestingly, one finds that /3 sharply 
increases when the temperature approaches the Curie 
temperature. This can be easily understood from the 
pronounced decreases of the spin relaxation time and 
the Zeeman splitting in this regime [see Figs.QJa) and 
[3]. By comparing the results with different values of 
the exchange coupling constant, one finds that (3 from 
J£x = — 1 eV is generally about one half of that ob- 
tained from Jpp — —0.5 eV because of the larger Zeeman 
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FIG. 5: (Color online) /3 and a as function of temperature 
with Jpp = -0.5 eV and -1.0 eV in sample A-C. In (b) and 
(d), the dots represent the experimental data from ferromag- 
netic resonance measurement for [001] (brown solid upper tri- 
angles), [110] (orange solid circles) , [100] (green open squares) 
and [1-10] (black open lower triangles) dc magnetic- field ori- 
entations fR.ef. l21T) . The brown solid square in (e) stands for 
the experimental result from domain-wall motion measure- 
ment (Ref.[H). 



splitting. Moreover, one notices that the nonmonotonic 
temperature dependence of the hole spin relaxation time 
in Fig.[3][b) is not reflected in /3 due to the influence of 
the Zeeman splitting. In all cases, the values of (3 can 
exceed one very near the Curie temperature. 

The results of the Gilbert damping coefficient from 
Eq. (fT3|) are shown as curves in Fig.[5][b), (d) and (f). 
The dots in these figures are the reported experimental 
data from the ferromagnetic resonance along different 
magnetic-field orientations. 21 Both the magnitude and 
the temperature dependence of our results agree well with 
the experimental data. From Fig.[2J one can conclude 
that the prefactor in Eq. (fT3)) . J/i/(iVMn|(S)|), is almost 
independent of temperature. Therefore, the temperature 
dependence of a mainly results from the nonadiabatic 
spin torque coefficient j3. Specifically, a is insensitive 
to the temperature in the low temperature regime and 
it gradually increases with increasing temperature due 
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to the increase of j3. Moreover, we predict that a be- 
gins to decrease with increasing temperature once (3 ex- 
ceeds one. This crossover lying at f3 « 1 can be expected 
from Eq. By comparing the results with different 

values of Jpp, one finds that the value of a is robust 
against the exchange coupling constant in the low tem- 
perature regime. In this regime, ^< 1 and one can sim- 
plify the expression of the Gilbert damping coefficient as 

a Ri jv M h s d (r app) • Since the total hole spin polariza- 
tion is proportional to the Zeeman splitting (see Fig. [2]) 
and t s is only weakly dependent on the Zeeman split- 
ting (see Fig. [3]) in this regime, the increase of Jg£ does 
not show significant effect on a. However, at high tem- 
perature, the scenario is quite different. For example, 
one has the maximum of the Gilbert damping coefficient 
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FIG. 6: (Color online) /3 and a as function of temperature by 
taking 7v M n = 5x 10 20 cm' 3 with J pp = -0.5 eV and -1.0 eV 
in sample D. The dots are from ferromagnetic resonance mea- 
surement (Ref.ll8l) for [001] (brown solid upper triangles) and 
[110] (orange solid circles) dc magnetic-field orientations. 

Since the effective Mn concentration of sample D is 
unavailable as mentioned above, we here take N-^n — 
5 x 10 20 cm~ 3 . The results are plotted in Fig.Hl It is 
seen that the Gilbert damping coefficients from our cal- 
culation with Jpp = — 1 eV agree with the experiment 
very well. As reported, the damping coefficient in this 
sample is much larger (~ 0.1) before annealing^ The 
large Gilbert damping coefficient in the as-grown sample 
may result from the direct spin-flip scattering between 
the holes and the random Mn spins, existing in low qual- 
ity samples. In the presence of this additional spin-flip 
channel, the hole spin relaxation time becomes shorter 
and results in an enhancement of a and j3 (for /3 < 1). 
Moreover, in the low temperature regime, a decrease of 
the Gilbert damping coefficient was observed by increas- 
ing temper aturejis which is absent in our results. This 
may originate from the complicated localization or cor- 
relation effects in such a disordered situation. The quan- 
titatively microscopic study in this case is beyond the 
scope of the present work. 

In addition, one notices that j3 in Ref.[U was deter- 
mined to be around 0.01, which is one order of magni- 
tude smaller than our result. The reason is because of 
the incorrect parameter used in that work, as pointed out 
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FIG. 7: (Color online) Spin stiffness (vertical spin stiffness) 
coefficient as function of temperature is plotted as curves 
with (without) symbols. The calculation is carried out with 
JS — —0.5 eV in sample E. The effective mass is taken to be 
1.0 (0.5)mo as labeled in the figure. The brown solid (from 
the period of the domains) and open (from the hysteresis cy- 
cle) squares are the experimental data of spin stiffness from 
Ref.Ea 



D. Spin stiffness and vertical spin stiffness 

In this subsection, we calculate the spin stiffness and 
vertical spin stiffness coefficients according to our previ- 
ous derivation 10 



and 



A ss = N h /[4m*{l+(3 2 )} 



A v ss = N h f3/[4m*(l + f3 2 )}. 



(15) 



(16) 



Since the effective mass to* is a rough description for the 
anisotropic valence bands in the presence of a large Zee- 
man splitting, it is difficult to obtain the accurate value 
of the stiffness coefficients from these formulas. Nev- 
ertheless, one can still estimate these coefficients with 
the effective mass taken as a parameter. The results are 
plotted in Fig. [7] By fitting the DOS of the occupied hole 
states, we find to* w too, which is consistent with the pre- 
vious work^ The spin stiffness and vertical spin stiffness 
coefficients with Nh = 3 x 10 20 cm~ 3 (1 x 10 20 cm~ 3 ) are 
plotted as the red solid (blue dashed) curves with and 
without symbols, respectively. The sudden decrease of 
A ss originates from the increase of /3 in the vicinity of 
the Curie temperature (see Fig.[S]). Our results are com- 
parable with the previous theoretical work from 6-band 
model. 2 - As a comparison, we take m* = 0.5mo, which 
is widely used to describe the heavy hole in the low en- 
ergy regime in the absence of the Zeeman splitting. 51 
The spin stiffness becomes two times larger. Moreover, 
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A^ s is found to present a nonmonotonic behavior in the 
temperature dependence as predicted by Eq. (|16l) . 

In Fig. we also plot the experimental data of the 
spin stiffness coefficient from Ref.|45|. It is seen that these 
values of A ss are comparable with our results and show 
a decrease as the temperature increases. However, one 
notices that the experimental data is more sensitive to 
the temperature especially for those determined from the 
domain period in the low temperature regime. This may 
originate from the strong anisotropic interband mixing 
and inhomogeneity in the real material. 

In Ref.hji we have shown that the vertical spin stiff- 
ness can lead to the magnetization rotated around the 
easy axis within the domain wall structure by Aip = 
(-y/ 1 + /3 2 — l)//3 in the absence of the demagnetization 
field. For /3 = 1, Aip « 0.13tt, while Aip = -> for 
^<1. As illustrated above, f3 is always larger than 0.1. 
Therefore, the vertical spin stiffness can present observ- 
able modification of the domain wall structure in GaM- 
nAs system^ 

IV. SUMMARY 

In summary, we theoretically investigate the tempera- 
ture dependence of the LLG coefficients in ferromagnetic 
GaMnAs, based on the microscopic calculation of the 
hole spin relaxation time. In our calculation, we employ 
the Zener model with the band structure carried out by 
diagonalizing the 8x8 Kane Hamiltonian together with 
the Zeeman energy due to the sp-d exchange interaction. 
We find that the hole spin relaxation time can present 
different temperature dependences, depending on the ef- 



fective Mn concentration, hole density and exchange cou- 
pling constant. In the case with high Mn concentra- 
tion and large exchange coupling constant, the hole spin 
relaxation time can be nonmonotonically dependent on 
temperature, resulting from the different interband spin 
mixings in the large and small Zeeman splitting regimes. 
These features are proposed to be for the estimation of 
the exchange coupling constant or itinerant hole density. 
By substituting the hole relaxation time, we calculate the 
temperature dependence of the Gilbert damping, nona- 
diabatic spin torque, spin stiffness, and vertical spin stiff- 
ness coefficients. We obtain the nonadiabatic spin torque 
coefficient around 0.1 ~ 0.3 at low temperature, which 
is consistent with the experiment. As the temperature 
increases, this coefficient shows a monotonic increase. In 
the low temperature regime, the Gilbert damping co- 
efficient increases with temperature, which shows good 
agreement with the experiments. We predict that the 
Gilbert damping coefficient can decrease with increasing 
temperature once the nonadiabatic spin torque coefficient 
exceed one in the vicinity of the Curie temperature. We 
also find that the spin stiffness decreases with increasing 
temperature and the vertical spin stiffness can present 
a nonmonotonic temperature dependence, similar to the 
Gilbert damping. 
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